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We study the motion of a classical particle interacting with one, two, and finally an infinite chain of ID square wells 
with oscillating depth. For a single well we find complicated scattering behavior even though there is no topological 
chaos due to the absence of hyperbolic periodic orbits. In contrast, for two coupled square wells there is chaotic 
scattering. The infinite oscillating chain yields the generic transition to chaos, with diffusion in energy and in space, 
as the separation between wells is increased. We briefly discuss the relevance of our results to solid state physics. 

PACS: 05.45. +b, 03.20 



1 Introduction 

Most of our knowledge about molecular, atomic, and ele- 
mentary particle physics comes from the analysis of scat- 
tering experiments. The potential fields producing the 
scattering may be the result of one or many particles and 
may be autonomous or time-dependent. Time-dependent 
scattering is particularly important nowadays in investi- 
gations of transport properties in microscopic and meso- 
scopic physics, where systems are usually modeled by 
simple Hamiltonians ( see, e.g., ||]). In this type of 
systems the description is generally performed using the 
Schrodinger equation since the devices under considera- 
tion are assumed to operate predominantly in the pure 
quantum limit. However, in order to get the complete 
picture of motion, it is necessary to study also their clas- 
sical motion, and the analysis of the quantum-classical 
correspondence is especially interesting when the classical 
motion becomes chaotic under certain conditions 

The study of the manifestations of chaos in the trans- 



port properties of mcsoscopic sytems has yielded impor- 
tant results. For example, measurements of magnetoresis- 
tance in ballistic microstructures in the shape of a chaotic 
billiard and circle showed clear distinctions in the power 
spectrum of universal conductance fluctuations and in the 
line shape of the weak localization peak In ||^ 

a criterion was proposed to distinguish between regular 
and chaotic dynamics by measuring classical resistivity in 
a type of mesoscopic channels. These results concern time 
independent phenomena. Here we report preliminary re- 
sults in the study of the classical dynamics of a type of 
simple time-dependent systems whose quantum counter- 
parts are models of microscopic or mesoscopic systems. 
Specifically, wc consider the scattering of classical particles 
interacting with one and two ID square wells with peri- 
odically oscillating depth. The square wells may represent 
the conduction band defined by a Al^Gai^xAs — GaAs 
heterostructurc or a quantum dot (see, e.g.,[0, ||ri|) and 
the oscillating depth may represent the electron-phonon 
interaction |12|] or the presence of a monochromatic elec- 
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tromagnetic or acoustic field , [|0[ ' llll ) • ^^^^ ^^^^ 
consider the dynamics of a classical particle interacting 
with a chain of oscillating wells. Since the chain is infinite 
and periodic the system is no longer an open system; the 
dynamics is recurrent in phase space. This system is the 
classical counterpart of a particle in a vibrating Kronig- 
Penney potential. The static Kronig-Penney potential in 
its quantum version has been fundamental in the under- 
standing of solid state properties (see, e.g., while its 
classical version is trivial. However, as will be shown be- 
low, when the chain vibrates the classical motion becomes 
chaotic. 

The organization of the paper is as follows. In sections 
n and HI, respectively, we study the scattering of particles 
off one and two oscillating wells. In Section IV, we consider 
the oscillating chain and in Section V we draw conclusions. 



2 A single oscillating well 

We consider classical particles interacting with the local- 
ized time-dependent potential. 



Vix,t) = 



Vo 

D cos((jjt) 



a; < 0, a; > a 
< a- < a, 



where a is the width of the well of height Vb, uj is the 
frequency of oscillation of the bottom of the well, and D 
is its amplitude (See Fig. |l|). We shall only consider the 
case D <Vo. Assume a particle approaches the well from 
the left, as indicated in Fig. |], with kinetic energy K^. 
At time t = tg the particle crosses the boundary a; = 
and the bottom of the well has a height Dcos{Loto). At 
that instant its kinetic energy changes abruptly to K* ~ 
Eq — DcosIujIq), where Eq = Kq + Vq. Since there is no 
gradient of x within the well, the "internal" kinetic energy 
K* is conserved while moving inside the well. After a time 
t = t* = a/p*, with p* = ^/l2K*) the particle reaches the 
right hand side of the well (a; = a), its total energy being 
El ~ K* + Dcos[LL>{to +t*)]. If El > Vo the particle exits 
to the right with energy 



where 



Ei= Ea + D{cos(j)i - cosM, (2) 



(f>i = 4>o + A(j)o, 00 = i^<o and (3) 



A0O = ujt* u}a/y^2{Eo - Dcos(j)o) (4) 

It is convenient to obtain dimensionless energies by re- 
scaling with respect to the height of the well Vq: 



ei ^ eo + rii{cos(j)i - cos(j)o) 
<Pi = 00 +V2/\/f^Q -riicos(j)o, 



(5) 
(6) 



_ 3_ 
-Dcos(mt)- 



Figure 1: Schematic of the ID oscillating well. Vb is the 
static depth and Kq is the kinetic energy of the incoming 
particle. 

where e = E/Vq, 



and 



772 = 2Truj /ujc, 
LUc = 2tt/Tc, and 
= a/^{2Vo). 



(7) 



(8) 



(1) Note that Tc is the time for a particle with total enery 
Vb to cross the well in the absence of oscillations. 

Equation (5) makes sense as long as ei > 1. If ei < 1, 
the particle reflects at a: = a without changing its internal 
kinetic energy K* . The particle then will reach a: = with 
energy and phase given by 



£2 = 

02 = 



£0 



77l(cOS02 - 

- A0O = 00 



:os0o) 
2A0O, 



where 



A0O = mN^a - Vicos(l)o 



(9) 



(10) 



Again, if £2 < 1 the particle will reflect towards x — a 
and so on, until for some number N of collisions with the 
boundaries, the energy cat becomes greater or equal to 1. 
The critical value eat = 1 implies that the particle exits 
the well but with zero kinetic energy. Thus, the scattering 
map is given by 



> 1 



eN = eo + rii{coscj)N ~ coscj)^) 
0JV = 00 + Nri2/ - Vicos(j>o (11) 

This map is deceptively simple. Clearly, if eo > 1 + 2ryi 
then ei > 1 and the particle will certainly be transmit- 
ted with phase and energy given by Eq.(ll) with = 1. 
However, if eo < 1 + 2771, then ei can be less than one, 
depending on the initial phase, and one cannot simply 
substitute into ejy to determine the fate of the scat- 
tered particle, i.e., whether it was transmitted or reflected 
and with what energy. It is necessary to calculate 0„ and 
e„ sequentially for n ~ 1,2,3 ... N to discover after how 
many collisions N the condition eAr > 1 is satisfied. 
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In order to get a panorama of the type of scattering 
asymptotes that can be obtained for typical values of 
rji and 772, it is useful to construct a hierarchy level 
plot where each level corresponds to a given number 
of collisions and the domain is the set of all incoming 
asymptotes JlSf . Such a hierarchy plot is shown in 
Fig. 2, for the parameters 771 = 0.5 and 772 = 16.45. 
This plot is typical of small to moderate frequencies of 
oscillations ( here 772 = 16.45, i.e., w ~ 2.6ujc). (For 
higher frequencies of oscillation of the well, the hierarchy 
plots show increasingly complicated structures). The 
boundaries between the level 1 and any other level N 
are given by the set of initial conditions (eo,0o) that 
yield outgoing energies eAr = 1. One may be tempted 
to obtain this set of (eo,^o) by inverting map (11) and 
using ejv = 1, but knowledge of the final phase (pN is 
lacking. However, the following procedure enables us to 
obtain this set. First we apply map (11) to the conditions 
eg = 1 and any (/jo in — tt < (^0 ^ This produces the 
outgoing values (e„,0„) for some integer n. Then it is 
easy to show, considering the symmetry of the cosine 
function, that forward iteration of (e„,— (/>„) yields the 
desired result e = 1. Hence, the initial conditions on the 
boundaries, or separatrices, between hierarchy levels are 
given by the set (cn,— f/'n) thus defined. Figure 3 shows 
the result of this procedure for a set of 6,000 values of 
initial phases in the range — tt < (I)q < tt. 
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Figure 2: Hierarchy level plot for scattering off one os- 
cillating square well. Each colored area represents the set 
of initial conditions that undergo a number of colli- 
sions (transversals) within the well before exiting. Light 
blue corresponds to TV greater or equal to 6. The white 
background is the set of initial conditions that transmit 
without any reflection. The parameters are: 771 = 0.5 and 
?72 = 16.45. 

With the help of hierarchy plots like Fig. 2 one can 
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Figure 3: Boundaries between hierarchy levels for the sin- 
gle oscillating well. Same parameters as in Fig. |[ See 
text for an explanation on the construction of this plot. 
771 = 0.5 and 7/2 = 16.45. 

easily determine for a given initial condition the outgo- 
ing energy and phase values, without having to iterate 
the map (11). One simply needs to locate the incoming 
asymptotes within its hierarchy level N and calculate 
Eq.(ll) for the known value N. That is, knowledge of the 
initial conditions enables us not only to determine, but to 
predict, the scattering values for N, e and (j). However, 
for sufBciently low values of energy, the structure of the 
hierarchy levels in general becomes more intricate, requir- 
ing higher precision in the specification of the incoming 
asymptotes. Figs. 4a and 4b show, respectively, a plot 
of the number of collisions (transversals) and a plot of the 
final energy as a function of (jjo for a fixed low incoming en- 
ergy. The great sensitivity to initial conditions is further 
demonstrated by Figs. 5a-5d. Figure 5a is an example 
of how the number of collisions can gradually increase ( 
here in steps of 2) as (j)o is varied over smaller and smaller 
intervals and then suddenly drops to = 1, direct trans- 
mission. Fig. 5b shows the corresponding behavior of the 
final energy. Fig. 5c is an example of a typical situation 
where changes in the initial condition yield unpredictable, 
seemingly random, changes in the number of collisions and 
Fig. 5d is the corresponding behavior of the final energy. 

The situation seems even more complex when the am- 
plitude of oscillation is as large as the height of the static 
well, (771 = 1). Figures 6 shows the rich structure of the re- 
sponse functions N{(I)q) and the outgoing energy e(0o)i due 
to the extreme sensitivity to initial conditions. Notice that 
this structure is nevertheless non-fractal. Numerically, the 
non-fractal structure is evident because of the absence of 
finer structures after increasing the resolution ( within nu- 
merical precision) enough times (see Fig. 7). Analytically, 
the fractal structure of the response ( or scattering) func- 
tions is ruled out by realizing that there arc no hyperbolic 
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Figure 4: Scattering functions for the single oscillating 
well ( 771 = 0.5,772 = 16.45). a) Number of collisions as a 
function of initial phase and b) final energy as a function 
of initial phase. The initial energy for all phases is eg = 
1.025. 



periodic points of the map. To see this consider the peri- 
odic orbits of period 1 . The condition for period one orbits 
(01 ^ 00, ei = eo) gives p* = m G Z. Within the 



well, p*^/2 + r]icos{4>) < 1, or, cos{4>) < 



'V2 
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Hence, 



one can find an integer m for a given set of parameters 
{a,u}, and 771) that satisfies the last inequality. This im- 
plies that there is a whole interval of the angle that 
is periodic of period one. Furthermore, since the period 
one points come in one-parameter families, i.e., they are 
not isolated, these points cannot be hyperbolic nor ellip- 
tic. Therefore, they must be parabolic ( see, e.g., [^). A 
similar argument rules out the existence of higher periodic 
hyperbolic points. Topological chaos appears in systems 
having transverse homoclinic or hcteroclinic connections 
(as this implies the existence of a Smalc horseshoe) which 
in turn can exist if there arc hyperbolic points [|l7[ . Con- 




Figure 5: Scattering functions for the single oscillating 
well (771 = 0.5,772 = 16.45) in two small ranges of initial 
phase. The ranges shown here correspond to magnifica- 
tions of Figures 4a and 4b. Note the two different types 
of behavior, namely gradual change ((a) and (b)) and ap- 
parently random change ((c) and (d)) in the number of 
collisions as a function of initial phase. 
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sequently. there cannot be topological chaos in the case 
of a single well. Finally, if there is no topological chaos, 
there is no chaotic scattering. 
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Figure 6: Scattering functions for the single well with rji ~ 
1.0, ?72 = 16.45. a (b) is the number of collisions ( final 
energy) as a function of initial phase. The initial energy 
for all phases is eg = 1.025. 



3 Two oscillating wells 

Now we analyze the interaction of a particle with a po- 
tential formed by two neighboring wells, as shown in Fig. 
8, oscillating in synchronization. The separation between 
the wells & is a new parameter that it is convenient to 
make dimcnsionless in the same way as was done for the 
well width a, namely. 



(12) 



Assuming the particle approaches the two wells from 
the left and arrives at a; = with energy eo and phase 
it will exit the first well with energy and phase given by 
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Figure 7: Magnification of the scattering function of Fig. 
6. 
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Figure 8: Schematic of the oscillating double well poten- 
tial. 



= eo + rii{cos(j>'„^^ - cos$o) 
= $o + ?7iiA$o, 



(13) 



where mi stands for the number of transversals within 
well number one. The reason for using the new symbols, 
primed phase 0^^^ and capitalized phase f&Oi will become 
clear shortly. If toi is even, the particle reflects from the 
first well and leaves the scattering region with energy emi 
but if mi is odd, the particle goes on to the second well. 
When it arrives ai x = a + b the phase is 



where 



S(t)mi = ml \/^nii - 1- 



(14) 



(15) 



is the phase increase while the particle travels from x — a 
to X = a + h. Note that the primed phase labels the value 
of the phase right when the particle exits a well and the 
capitalized phase labels the phase as it enters the well. We 
shall use this convention throughout the rest of the paper. 
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After transversing the second well m2 times the particle 
exits with energy and phase given by 

C = «'™i+m2A$„,. (16) 

If 7712 is odd the particle exits the second well at x = 
a+b+a, leaving the interaction region, but if m2 is even the 
particle exits at a; = a + 5, returning toward the first well, 
and so on. Let rrii stand for the number of transversals 
within the well number i, where i odd (even) refers to well 
number one ( two). We can then write 

em„ = em,._i + ?7i (cos0'„^^ - cos4>„^_ J 

?^'m„ = $m„_i +m„A$™„_,, (17) 

where 

A<I>„,^ = 772 /Vem. - ?ncos$„, (18) 

and nii is even for all i except i = 1 and i ~ n { The 
case mi= even is given by Eq. (13), where the parti- 
cle did not interact with the second well), n labels the 
last well the particle crossed before leaving the interact- 
ing region and so is also the total number of interactions 
with the wells. On the other hand, the total number N 
of transversals is the sum of all irii from i = 1 to « = n. 
Figure 9 is the hierarchy plot computed for 773 = 5 and 
with 771 and 772 the same as for figure 2. Note that for two 
wells, a total number of A'' = 2 transversals can be pro- 
duced in two ways: either, two transversals in the first well 
with no interaction with the second well, or one transver- 
sal in the first well and one with the second well. Figure 
9 associates the color blue to the latter possibility while 
the white area corresponds to the first possibility. Com- 
parison with Fig. 2 shows that the smooth boundaries 
between various hierarchy levels now appear fragmented 
and the scattering functions are expected to show some 
evidence of this fragmentation. Figures 10a and 10b plot 
the energy and the number of collisions, respectively, as 
a function of the initial phase for the same parameters as 
those of Fig. 9 and for an incoming energy eq = 1.025. At 
first sight the plots N{(j)o) do not appear qualitatively dif- 
ferent from those corresponding to a single well (Figs. 4 
and 6) but a closer look reveals that as the resolution is 
increased there appears more and more structure, suggest- 
ing its fractality ( see Fig. 11). That such fractality exists 
can be demonstrated by considering the dynamics of par- 
ticles initially placed between the wells. Let us first define 
the Poincare surface to be the midpoint between the two 
wells and positive momentum and choose the coordinates 
to be the energy and phase. The energy at the Poincare 
surface is given by the first of Eqs. (17) but the phase 
is shifted by half of <5(/)m„ to correspond to the midpoint 
between the wells. Thus, the Poincare map is 




Figure 9: Hierarchy level plot for the oscillating double 
well. There are two ways of producing A^ = 2 as mentioned 
in the text. The white area corresponds to direct trans- 
mission (one transversal in the first well and one in the 
second one). The blue area corresponds to two transver- 
sals in the first well and none in the second one. In general, 
there are 

2N/2 

ways. Here all possibilities are plotted for 
A^ equal or greater than 4. The parameters are jyi = 0.5, 
'72 = 16.45, and 773 = 5. 

We now place a set of initial conditions on a grid on this 
surface and let them evolve forward in time. Some initial 
conditions will leak after a finite number of interactions 
with one or both wells and some will remain trapped for- 
ever. Fig. 12a is a plot for those orbits that return to the 
Poincare surface at least once and up to 1000 times. In 
Fig. 12b we plot only those orbits that returned to the 
surface at least 100 times and up to 1000 times. That 
is, transient orbits that crossed the Poincare surface less 
than 100 times were not plotted in Fig. 12b. Finally, in 
Fig. 13 we plot orbits that cross this Poincare section 
at least 4000 times and up to 5000 times. This Poincare 
plot now represents the dynamics of the particles that are 
trapped forever (bounded orbits) and shows the appear- 
ance of an island of stability associated with a fixed point 
of period one and the surrounding elliptic islands of higher 
period ( here period 8). This is just the common structure 
predicted by the Poincarc-Birkhoff theorem ||l8| with al- 
ternating elliptic and hyperbolic fixed points. The hyper- 
bolic points and their associated homoclinic/heteroclinic 
tangle are not shown in this figure. The outermost island 
chains are not screened behind a KAM curve, therefore 
the scattering flow can reach the homoclinic tangle. Cor- 
respondingly, the tendrils of the stable and unstable man- 
ifolds of this tangle reach out into the asymptotic region, 
thus producing the fractal structure of the scattering func- 
tions ijl^. Hence, we have a numerically based proof that 
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Figure 10: Scattering functions for the oscillating double 
well, a) Number of collisions as a function of initial phase 
and b) final energy as a function of initial phase. The 
initial energy for all phases is eg = 1.025. rji = 0.5, r/2 = 
16.45, and 773 = 5.0 
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Figure 11: Magnification of Figure 10b. 

points have been labeled (e„, (?!)„) , n — 1,2,3... instead of 
{^m^,4>m„), n = 1,2,3..., because these plots do not con- 
tain the information on how many transversals m„ occur 
for each interaction n. 

For the analysis of these plots it is convenient to write 
the energy e„j^ of eq. (19) in terms of the initial condi- 
tions: 

em„ = + rii{cos(t>[^^ ~ cos<^q) + 

n-l 



=1 



(20) 



Another formula that will be useful below is obtained 
by substituting (j)'^^ = (j)m„ - ^S(j)m„ and = 0m„ -I- 
i(50„i„ in (20) and expanding the cosine terms: 



a particle scattering off this double well oscillating poten- 
tial gives rise to true chaotic scattering, unlike the case of 
a single well. 

4 Oscillating chain 

The infinite oscillating chain is formed by the infinite rep- 
etition of the double well arrangement and hence (19) is 
the Poincare map describing the motion of the particle, 
except that m.i has a different meaning now. Namely, 
m.i now stands for transversals during the ith inter- 
action of the particle with a well. Of course, here 

can be odd or even for any i, i ~ 1,2,3 Since the 

chain is periodic the Poincare plots are modulo a + b. Fig- 
ures 14-17 show the Poincare plots for increasing values 
of ?73 (r/i — 0.5; r/2 ~ 16.45). In these figures the phase 



em„ = £0 + ?/l[cOs(0m„ - -^Hni„) - COS$o)] + 
n — 1 ^ 

+ 2ryi ^ sin(j)m, sin-6(j)mi (21) 



+ ^'50„,„_, (22) 

Some features of the Poincare plots for small 7/3 can 
be understood qualitatively by first analyzing the un- 
perturbed case ( 773 = 0). If 773 = then 5(t)mi = 0, 
= 4'mi — 4>'m- f^i' and (21) and (22) reduce to 

em„ = Co + ??i(cos0„„ - cos^o), (23) 
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Figurc 12: Transient and bound orbits for the oscillating 
double well, a) is the Poincare plot for those orbits that 
return to the surface once and up to 1000 times; b) same 
as a) but plots only those that crossed the surface at least 
100 times and up to 1000 times. Here 772 = 6.28, 773 — 0.1, 
and r\\ = 0.5 

n 

= 0o + A</)o^m,. (24) 

1=1 

To obtain the last equation we used the fact that 
= A(/)o for all i since ~ rjicoscjimi is the (adi- 
mensional) internal kinetic energy, which is conserved for 
773 = ( see discussion below Eq. (1)). Eqs. (23) and (24) 
show that for 773 = all orbits are described by invariant 
ID curves, the cosine curves em„(eo,0o)- Recall that for 
a given interaction n with a well e,„^ > 1 but e,„/^ < 1 for 
mjj < 77i„. It is clear from (23) that if eq > 1 + 2771 then 
> 1 and 77i„ — 1 for all n and hence the motion is 
of the translational type; the particle ncvcrs changes di- 
rection. If eo = 1 + 2771 then the right hand side of (23) 
can be as small as one for some n, provided 00 = and 



1.17 
1.16 




Figure 13: Poincare plot for the oscillating double well. 
Here the minimum number of crossings can be as large 
as 5000 for some orbits. Orbits that left the interactions 
region before 4000 crossings are not plotted. Same param- 
eters as Fig. 12. 

0m„ = kiT, fc, odd. This implies that 7iA0o — 77772 + Vi 
must be exactly equal to kn, but for any 771 in the range 
(0, 1) and any 772 > 0, there are some integers n and k that 
approximate arbitrarily close this equality. Since e = 1 is 
the energy threshold, above which the particle does not 
reflect, then it follows that the curve e = 1 -I- 7^1 -I- rjicoscf) 
defines the separatrix curve, separating translational or- 
bits from librational orbits ( curves below the separatrix). 
For initial conditions placed below the separatrix, m„ is 
in general larger than one. Since all curves, translational 
or librational, are invariant curves, there is no possibility 
for a translational orbit to become a librational one, and 
viceversa. The three top orbits of Fig. 14a are transla- 
tional and the rest are librational. 

Now consider the 773 > case. Notice that the last term 
of (21) produces, as a function of (j>m„ — 5'50?ri„, fluctua- 
tions from the invariant cosine curve. These fluctuations 
can be small or large, depending on the value of 773 and 
whether the terms in the sum add constructively or de- 
structively. Consider the case 7/3 << 1. If in addition 
^(^mi << 1 for all i then the summation term can be 
approximated by 771 ^ sin(j)miS4>rni- The summed terms 
would yield an appreciable deviation from the cosine curve 
if sin(j)mi were to oscillate slowly even if each term in the 
sum is small. This could happen for low energies ( c.f., 
(22) and (18)) and the orbit in the Poincare surface will 
appear as scattered points in the neighborhood of the co- 
sine curve. Similarly, if e„i. is close to one, the last term 
of (19) may become the dominant term. Thus even for a 
small 773, the perturbation can be strong for initial condi- 
tions in the neighborhood of the unperturbed separatrix or 
for initial conditions on low-lying unperturbed librational 
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orbits, see Fig. 14b. This is the concrete mechanism in 
this system for the appearance of chaotic orbits as 7/3 be- 
comes sHghtly larger than zero ( Figs. 14b and 15). 

Further analysis on the development of chaos as 773 
is increased will be reported elswhere. Here we shall 
only remark on some features observed in our numerical 
results ( Figs-15-17). It is interesting to note, for example, 
that the "boomerang" like stable island immersed in the 
chaotic sea of Fig. 15b ( at e„ « 1.2) corresponds to 
the period one stable fixed point of the double well. In 
fact, the Poincare plot for orbits trapped between two 
wells (see Fig. 13) fits exactly within the abovementioned 
stable island of the chain, and the heteroclinic tangle 
of the double well lies on the chaotic sea of the chain. 
The elliptic islands above the lowest one most likely 
correspond to period one stable fixed points produced by 
trajectories trapped between three, four, and more wells. 

Figures 15b to 17 show that as r/3 continues to increase, 



the chaotic area increases, allowing the energy to diffuse 
to higher and higher values. Thus if the energy is initially 
very small, after some time it may become very large as 
long as diffusion in energy is not limited by KAM curves. 

In addition to energy diffusion, there is also diffusion in 
configuration space since the particle is interacting chaot- 
ically with the wells. Figure 18a plots the well number vs. 
number of interaction for the same parameters as in Fig. 
17. Notice the random- walk- like behavior. That is, the 
particle, starting at well number one. seems to travel ran- 
domly forward and backward, returning repeatedly to the 
starting point. While this is happening, the energy of the 
particle oscillates wildly between one and approximately 
seven (sec Fig. 18b). Then after approximately 55,000 in- 
teractions ( see Fig. 18b), the particle travels only in the 
negative direction during several interactions. This corre- 
sponds to a jump in its energy and the particle diffuses 
to higher energy values ( sec Fig. 18b). This behavior 
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Figure 16: Poincarc plots for the oscillating chain of wells, 
a) 773 = 0.1 and b) 773 = 10. Same 771 and 772 as in Fig. 15. 
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Figure 17: Poincarc plot for the oscillating well with 773 
500. Same r/i and 772 as in Fig. 16. 



can be partially understood by examing the Poincarc plot 
of Fig. 17. Notice that for energies around seven units, 
there are some visible white areas indicating the existence 
of some large, very thin, islands of stability, where the 
orbits get stuck temporarily. For larger number of inter- 
actions the particle eventually returns to the original well, 
performing a random-like-walk and the energy fluctuates 
over the whole range delimited by the KAM curve found 
around e = 42. 
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Figure 18: Space and energy diffusion in the oscillating 
chain, a) is diffusion in space, represented by the well 
number and b) is diffusion in energy of the particle as 
it interacts with the various wells. The parameters are 
7/1 = 0.5, 7/3 = 16.45, and 773 = 500. 



5 Conclusions 

In this paper we have described the motion of a classical 
particle interacting with one, two, and finally an infinite 
chain of oscillating square wells. For the single well we 
found that although the scattering functions, such as final 
energy and number of collisions as a function of initial con- 
ditions, may yield complex behavior, extremely sensitive 
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to initial conditions ( hence unpredictable), the nonexis- 
tence of hyperbolic fixed points renders the scattering non- 
chaotic. On the other hand, when the potential region is 
formed by two oscillating wells, the scattering functions 
are truly fractal as expected since it was shown that there 
exists topological chaos. The dynamics of a particle in- 
teracting with a chain of oscillating wells was shown to 
yield the transition to chaos as the separation between 
the wells increases. Although there is chaotic diff'usion 
in energy this is limited by KAM curves and the diffu- 
sion in configuration space is random-like but with some 
peculiarities. The scattering off one and two oscillating 
square wells has been considered previously by some au- 
thors to model some microscopic and mesoscopic systems 
but their treatments have been completely quantum me- 
chanical. Here we treat the purely classical description 
as a starting point to understand their classical-quantum 
correspondence. The oscillating chain represents a clas- 
sical oscillating Kronig-Penney model, particularly rele- 
vant to solid state systems, and it may help in the un- 
derstanding of some quantum properties. For example, 
some localization of the wave function of an electron in 
a one-dimensional crystal or in a superlattice in the pres- 
ence of some monochromatic field ||2^ is expected to oc- 
cur since we have shown that the motion can be chaotic. 
This should manifest itself in the reduction of conduc- 
tivity. Thus, there is some similarity with the Anderson 
Model |2^. The dynamics of a particle in the oscillating 
chain is also similar to the Fermi accelerator model, with 
the additional nice feature that in the chain there is dif- 
fusion in space as well as in energy. Finally, as far as the 
energy is concerned the oscillating chain is equivalent to a 
single oscillating well inside a larger static infinite square 
well. An oscillating square barrier inside a larger static 
infinite square well has been recently considered in ||2^ . 
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